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ABSTRACT 

An important test for models of galaxy formation lies in the metallicity distribu¬ 
tion functions (MDFs) of spheroid stars and their globular clusters (GCs). We have 
compared the MDFs obtained from spectroscopy of the GCs and the star-by-star pho¬ 
tometry of the old halo red giants in the nearby elliptical galaxy NGC 5128, with the 
predictions of a ACDM semi-analytic galaxy formation model. We have selected model 
ellipticals comparable in luminosity and environment to NGC 5128, and reconstructed 
their MDFs by summing the total star formation occurring over all their progenitors. 
A direct comparison between models and data shows that the MDFs are qualitatively 
similar, both have stellar components which are predominantly metal-rich (~0.8^q), 
with a small fraction of metal-poor stars extending down to 0.002 zTq. The model 
MDFs show only small variations between systems, whether they constitute brightest 
cluster galaxies or low luminosity group ellipticals. Our comparison also reveals that 
these model MDFs harbour a greater fraction of stars at Z > Zq than the obser¬ 
vations, producing generally more metal-rich (by ~0.1 dex) MDFs. One possibility 
is that the outer-bulge observations are missing some of the highest-metallicity stars 
in this galaxy. We find good agreement between the model and observed GC MDFs, 
provided that the metal-poor GC formation is halted early (z ~5) in the model. Under 
this proviso, both the models and data are bimodal with peaks at O.IZq and Zq, and 
cover similar metallicity ranges. This broad agreement for the stars and GCs suggests 
that the bulk of the stellar population in NGC 5128 may have been built up in a 
hierarchical fashion, involving both quiescent and merger-induced star formation. The 
predicted existence of age structure amongst the metal-rich GCs needs to be tested 
against high-quality data for this galaxy. 

Key words: galaxies: interactions - galaxies: elliptical - galaxies: evolution - galaxies: 
individual: NGC 5128 - galaxies: star clusters 


The task of understanding the star formation histories of 
elliptical galaxies is one of the most challenging areas of 
galaxy formation, and has generated a wide and steadily 
growing literature. For giant elliptical (gE) galaxies, the for¬ 
mation history may be particularly complex because these 
systems are likely to involve some (if not all) of the major 
formation processes that are currently debated, including 


1 INTRODUCTION 


singular or multi-phase collapse (e.g. Silk 1977; Arimoto & 
Yoshii 1987; Forbes, Brodie, & Grillmair 1997), the hier¬ 
archical merging of gaseous fragments in the early universe 
(e.g. Kauffmann, White, & Guiderdoni 1993; Dubinski 1998; 
Pearce et al. 1999; Cole et al. 2000; Somerville, Primack, & 
Faber 2001), the major merging disc galaxies (e.g. Schweizer 
1987; Ashman & Zepf 1992; Barnes & Hernquist 1992; Sh- 
ioya & Bekki 1998; Hibbard & Yun 1999; Naab & Burkert 
2001; Bekki et al. 2002) and ongoing dissipationless accre¬ 
tion of satellites (e.g. Cote, Marzke, & West 1998; Cote, 
West, & Marzke 2002). 
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Theoretical analyses employing sophisticated N-body 
simulations which include chemical evolution (e.g. Shioya 
& Bekki 1998; Steinmetz & Navarro 1999; Koda, Sofue, & 
Wada 2000; Kawata 2001; Lia, Portinari, & Carraro 2002) 
and semi-analytic treatments (e.g. Cole et al. 1994, 2000; 
Kauffmann 1996; Somerville & Primack 1999; Thomas 1999) 
are rapidly adding to our ability to understand this early se¬ 
quence of events. 

However, observational constraints are badly needed to 
sort out the relative importance of these possibilities. In this 
paper, we explore the ability of one hierarchical merging 
code to reproduce the metallicity distributions of the re¬ 
cently obtained star-by-star halo photometry and GC sys¬ 
tem of the nearby elliptical galaxy NGC 5128. 

At high-redshifts, there is some observational evidence 
that the bulk of the stars in at least a subset of galaxies may 
be in place as early as 2 > 5 (e.g. Dunlop et al. 1996; Franx 
et al. 1997; Dey et al. 1998; Nolan et al. 2001; Waddington 
et al. 2002). These studies have been generally used to argue 
for single-collapse, high-redshift models of gE galaxy forma¬ 
tion. Other discussions suggest a somewhat longer timescale 
and later epochs for much of the star formation (e.g. Zepf 
1997; Kauffmann & Chariot 1998; Barger et al. 1999; Menan- 
tean, Abraham, & Ellis 2001 - but also see e.g. Jimenez et al. 
1999 for an different interpretation of the Zepf 1997 results). 

On the observational side, an opportunity to test these 
models is through the metallicity distribution function 
(MDF) of the oldest stars in ellipticals. Indirect estimates 
of the MDF can be made through measures of the spectral 
features of their integrated bulge light (e.g. Gonzalez 1993; 
Davies, Sadler, & Peletier 1993; Kuntschner & Davies 1999; 
Trager et al. 1998; Terlevich & Forbes 2002). However, two 
considerably more direct ways to study the MDF are now 
emerging as useful tools in the low redsliift regime: star-by¬ 
star photometry of the old-halo red giant stars; and metal- 
licity measurements of the old-halo globular clusters (GCs) 
within the galaxies. Despite observational biases and short¬ 
comings (see below), these kinds of direct object-by-object 
construction of the MDF provide a level of new information 
that cannot be obtained in any other way. 

Neither of these latter approaches is immune from the 
interpretive difficulties associated with all such studies: for 
example, metallicity measurements of either the old red- 
giant stars or GCs are uncertain to the degree that we know 
the particular mixture of ages present in the system, or even 
whether a well defined age-metallicity relation (AMR) ex¬ 
ists. And for the GCs, other issues arise such as the effi¬ 
ciency of cluster formation, which might in principle vary 
with metallicity or even the epoch of formation. 

Recently, Beasley et al. (2002) have applied a semi- 
analytic galaxy formation code (Cole et al. 2000) to synthe¬ 
sise model MDFs for globular cluster systems (GCSs), and 
have compared these with real GCSs in early-type galax¬ 
ies. Detailed discussion of the elements of the model can be 
found in these two papers. In the model, a large galaxy is 
progressively assembled from a population of highly gaseous 
proto-galactic fragments (in the Cole et al. model, these are 
modelled as rotationally supported discs, and in this study 
we use the terminology of proto-galactic discs (PGDs) used 
by Beasley et al. 2002). The most important feature for our 
purposes is that the formation process, during which time 
star formation goes on simultaneously with the hierarchical 


merging of these PGDs, gives rise to two main subpopula¬ 
tions of objects: 

(i) a population of stars and metal-poor GCs that are built 
within the PGDs before they merge into the body of the 
larger central galaxy (termed the “quiescent mode” of star 
formation); and 

(ii) a predominantly metal-rich population built in suc¬ 
cessive star-bursts generated during major mergers (“burst 
mode”). This burst mode can happen any time that a PGD 
(in the fiducial model defined as a minimum of 30% of the 
central galaxy mass) merges with the central galaxy, and so 
can continue until quite late epochs. 

The quiescent mode of star formation generates a dis¬ 
tribution over age and metallicity that is fairly smooth and 
continuous (star formation in this mode begins at 2 ~ 12 
in the fiducial model, and continues at low levels until late 
epochs). The burst mode also produces stars over a signifi¬ 
cant fraction of a Hubble time, although the onset of merg¬ 
ing begins after the first PGDs have formed from cooling 
gas. In this mode, the stellar populations formed are more 
metal-rich due to their formation from enriched gas, and are 
in the mean younger. Because these mergers are stochasti¬ 
cally driven, the age and metallicity distributions formed in 
bursts can differ widely from one galaxy to another. 

It is natural to associate the bimodal form of the GC 
MDF seen frequently for gE galaxies with the two subpop¬ 
ulations described above (e.g. Zepf & Ashman 1993; Forbes, 
Brodie, & Grillmair 1997; Larsen et al. 2001). Typically, 
more or less equal numbers of GCs fall within each mode, 
and the two modes have metallicity peaks consistently near 
[Fe/H] ~ —1.5 and ~ —0.5 roughly independent of galaxy 
luminosity. In this interpretation, the “blue” GCs are seen 
as belonging to the earlier, quiescent mode while the “red”, 
more metal-rich ones arise from the burst mode. 

However, Beasley et al. (2002) find that to predict the 
right numbers and metallicities of the GCs in each mode, 
it is necessary to impose two extra numerical factors: first, 
the formation times of the blue GCs formed in the quiescent 
mode must be truncated at an early epoch Ztrunc (other¬ 
wise they would continue forming to later and later times at 
higher and higher metallicities, and their MDF would then 
end up overlapping substantially with the red GCs). Second, 
the relative formation efficiency, defined as the mass ratio 
of GCs to all stars, e = Mgc/M *, needs to be different for 
the two modes (otherwise the total numbers of GCs in each 
mode would be widely different). For the well studied proto¬ 
type case of NGC 4472, Beasley et al. (2002) find best-fitting 
values Ztrunc — 5, cr = 0.007, and es = 0.002. As we show 
later (Section 6) the numerical value of cr is sensitive to a 
least one parameter in the semi-analytic model, and should 
not be taken as an ’absolute value’. The underlying physi¬ 
cal reasons for enforcing the truncation at 2 ~ 5 is not yet 
understood, but within the parameters of the semi-analytic 
model, it is necessary to impose this condition to maintain 
the observed numbers of GCs and the form of the MDF. 

An obvious way to develop these tests further is to com¬ 
pare the predictions of the semi-analytic model with the 
MDFs for both the GCs and field-halo stars simultaneously 
in the same galaxy. A new opportunity to do this has arisen 
with NGC 5128, the giant galaxy at the centre of the Cen- 
taurus group at d ~ 4 Mpc (Harris, Harris, & Poole 1999). 
Through several recent HST-based studies, its halo red-giant 
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Figure 1. ; Upper panel: Observed [m/H] distribution for ~ 17,000 field stars in NGC 5128 at a location 8 kpc SW of the galaxy centre, 
from Harris &; Harris (2002). This sample of stars is completely dominated by an old red-giant population and has relatively high mean 
metallicity, with very few stars more metal-poor than [m/H] ~ —0.8. Lower panel: Observed [m/H] distribution for GCs in the halo of 
NGC 5128, from the spectroscopic sample of Held et al. (2002). 


Table 1 . Some basic properties of NGC 5128 and our selected model galaxies. 


Quantity 

NGC 5128 

Galaxy #1 

#2 

#3 

#4 

#5 

#6 

#7 

#8 

M s 

-21.0“ 

-22.63 

-21.49 

-21.28 

-21.13 

-21.06 

-20.97 

-20.61 

-20.32 

<t( kms —1 ) 

114 6 

1460 

368 

286 

330 

244 

288 

247 

228 

m/c b 

3.9“ 

3.0 

3.2 

4.3 

5.5 

3.2 

6.6 

5.1 

4.0 

(B - V) 0 

0.88“ 

0.89 

0.92 

0.89 

0.91 

0.91 

0.91 

0.89 

0.89 

(V - K)o 

2.79 d f 

2.93 

3.02 

2.91 

2.97 

2.97 

3.00 

2.94 

2.93 

M m (x 1O 9 M 0 ) 

0.83T0.25 “ 

3.1 

2.1 

1.0 

1.3 

0.49 

1.5 

0.065 

1.3 

S N 

2.6±0.6 / 

5.8 

4.1 

3.8 

2.8 

3.7 

2.5 

2.0 

2.7 

p( KS)9 

1.0 

0.35 

0.63 

0.66 

0.80 

0.78 

0.63 

0.78 

0.35 


“ B-band absolute magnitude, integrated B — V colour (de Vaucouleurs et al. 1991), b Cen A group velocity dispersion (van den Bergh 
2000), “ B-band mass-to-light ratio (Hui et al. 1995), d f integrated V — K colour (Pahre 1999), calculated with unmatched apertures. “ 
Neutral gas mass (Richter, Sackett, & Sparke 1994). f GC Specific frequency (Harris 1991). 9 Kolmogorov-Smirnov probability. 


stars have been resolved and used to construct well defined 
MDFs (Soria et al. 1996; Harris, Harris, & Poole 1999; Mar- 
leau et al. 2000; Harris & Harris 2000, 2002). In addition, 
metallicity measurements based on spectroscopic indices are 
now available for ~ 40 of its halo GCs (Held et al. 2002). 
The time is thus ripe to combine this observational material 
to test the validity of these models, and see if a consistent 
formation history of this galaxy can be constructed. 

The paper is organised in the following way: In Sec¬ 
tion 2, we summarise the observational data used in this 
study. In Section 3, we briefly discuss the aspect of the semi- 
analytic model most relevant to this study, namely the star 


formation recipe employed. Next, in Section 4, we show ex¬ 
ample outputs of the MDFs and star formation histories of 
eight model elliptical galaxies and their associated GC sys¬ 
tems. In Section 5, we pursue a specific comparison between 
one arbitrarily selected galaxy and the observational mate¬ 
rial. In Section 6 we vary the parameter which controls the 
threshold for star-bursts in the model to test the robustness 
of the results from the fiducial semi-analytic model. In Sec¬ 
tion 7, we propose a possible age-metallicity relation for the 
NGC 5128 GCs to be tested against future data and finally 
in Section 8 we summarise our conclusions. 
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2 SUMMARY OF THE OBSERVATIONAL 

DATA 

For the stellar halo of NGC 5128, MDFs based on large 
samples of stars have been obtained through HST/WFPC2 
(V, I) photometry. A finely spaced grid of red-giant evo¬ 
lutionary tracks, calibrated against Milky Way GCs, is 
then superimposed on the colour-magnitude diagram for the 
NGC 5128 stars and used to define the stellar metallicity 
distribution (Harris, Harris, & Poole 1999; Harris & Harris 
2000, 2002). To date, MDFs are available for three different 
locations in the NGC 5128 halo. The first two are located 
at projected galactocentric distances of 21 and 31 kpc (Har¬ 
ris, Harris, & Poole 1999, Harris & Harris 2000) and so are 
likely to be representative of only the outer halo. The third 
of these is a field located at 8 kpc which therefore samples 
a combination of the halo and bulge (see the discussion of 
Harris & Harris 2002). For our present purposes, we will use 
this inner (halo plus bulge) sample to define a broadly based 
MDF representing a plausible average over the galaxy as a 
whole. 

In this paper we compare our model results for GCs 
with the spectroscopically derived metallicities of Held et al. 
(2002). These indices are transformed into [m/H] = log 
(Z/Zq) as described in Harris & Harris (2002), yielding ^40 
GCs drawn widely from all over the NGC 5128 halo. We note 
that the GC photometric metallicities from the (C — Ti) in¬ 
dices of Harris et al. (1992) and the (U — V) indices of Re- 
jkuba (2001) are consistent with the Held et al. (2002) GC 
metallicity distribution. All three of the above GC studies 
probe regions well beyond the central dust lane (which con¬ 
tains a small number of < 1 Gyr old GC candidates asso¬ 
ciated with a recent merger event, e.g. see Holland, Cote, & 
Hesser 1999). The [m/H] histograms for both the field-star 
and spectroscopic GC sample are shown in Figure 1. 

Although small-number statistics may still be affect¬ 
ing our estimates of the true proportions of metal-poor and 
metal-rich GCs (lower panel of Fig 1), the GC and field- 
star MDFs are clearly quite different. Roughly half the GCs 
are more metal-poor than [m/H]= 1, whereas remarkably 

small numbers of field stars populate this low-metallicity 
range. Another way to quantify this is to say that the specific 
frequency of GCs to field stars is a function of metallicity 
(Harris & Harris 2002). The stellar-halo MDF (upper panel 
of Fig. 1) has a single major peak at [m/H] ~ -0.4, though it 
is clearly not symmetric, having a long and thinly populated 
low-metallicity tail. These two very different MDFs - stel¬ 
lar and GC system - must be matched by the semi-analytic 
model simultaneously. 


3 STAR FORMATION IN THE MODEL 

The semi-analytic model we test in this study is the fidu¬ 
cial version of GALFORM, as described in detail by Cole 
et al. (2000)(where h=0.7, with Ho = 100 h kms -1 Mpc -1 , 
flo=0.3, I2b=0.7, Ao=0.7). Its application to the modelling 
of the GC systems of elliptical galaxies is presented in 
Beasley et al. (2002). 

This fiducial version of the Cole et al. model was devel¬ 
oped to best match local galaxy properties, in particular the 
local R-band luminosity function. As such, the model should 


possess genuinely predictive powers. Therefore, for the ma¬ 
jority of this study, we retain the identical parameter values 
given in Cole et al. (2000). However, in Section 6 we vary 
the parameter which controls the burst threshold to test the 
robustness of the form of the MDF in the model. 

Since we are examining the stellar metallicity distribu¬ 
tions of the model ellipticals here, it is the details of the 
model star formation that we are primarily interested in 
here. We briefly discuss the salient points, although for a 
full description see Cole et al. (2000). In the model, the 
PGDs are modelled as rotationally supported, exponential 
gas discs. This is for convenience only, whereas in reality 
they may be expected to form a continuum of objects from 
flattened discs, to more clumpy systems dominated by ran¬ 
dom motions. To avoid possible confusion with spiral discs, 
we continue to refer to these objects as “PGDs”. 

Star formation in the model proceeds in two modes, 
quiescently in PGDs and during the merging of these PGDs. 
The instantaneous star formation rate (ip) occurring in these 
PGDs (the quiescent mode) is proportional to the mass of 
cold gas available (M co ia) 

, Afcold / -. 

1P= - (1) 

T* 

where r* is the star formation timescale, which is defined as 

r* = e/VpGD (Vpgd/200 kms _1 )“* (2) 

here, t P gd = cpgd/Vpgd, where Vpgd and ?- P gd are the 

circular velocity and half-mass radius of the PGDs respec¬ 
tively, whilst ejT 1 and a* are dimensionless parameters in 
the model. 

Star formation also occurs during major mergers, con¬ 
trolled by the parameter / e m p , which is the mass ratio of 
the merger progenitors. In the fiducial model of Cole et al. 
(2000), f e nip = 0.3, corresponding to a minimum mass frac¬ 
tion of 30% required to induce a major merger. During these 
major mergers, star formation proceeds similarly to that oc¬ 
curring quiescently in the PGDs, but is based upon the prop¬ 
erties of the central spheroid (i.e. its dynamical time) rather 
than the properties of the merging PGDs. In the fiducial 
model used here a* =0, which as discussed in Cole et al. 
(2000), yields a r* oc tpgd scaling in equation 2, as sug¬ 
gested by results such as those of Kennicutt (1998). 

In lockstep with this star formation, stellar winds and 
SNe are assumed to ‘feedback’ energy and mass into the 
IGM, re-heating cooled gas and depositing metals into the 
hot gas reservoir of the halo. The level of mass deposition 
from these PGDs is governed by the star formation rate and 
the feedback efficiency (/?) 

dMeject = flip d t (3) 

with this feedback efficiency defined in terms of the PGD 
circular velocity 

(3=(V PGD/I4ot) _ “ hot (4) 

where ahot is a dimensionless parameter, and Vhot is a pa¬ 
rameter with units of kms -1 . 

In the above recipe, up to ~ 30% of the total stellar mass 
of the model galaxies is composed of burst-formed stars, the 
rest originate from the quiescent mode. This, however, is not 
the only star formation prescription used in semi-analytic 
modelling. Somerville & Primack (1999) present a discussion 
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Figure 2. Sample metallicity distribution functions for the stars in eight semi-analytic model galaxies. In each panel, the B-band 
magnitude (Mg) of the galaxy and the velocity dispersion (cr) of the halo in which it resides is labelled at upper left. The numbers in 
parenthesis in the upper identify the number of the model galaxy listed in Table 1. The light shaded part of each histogram shows the 
metallicity distribution for stars formed in the quiescent mode, while the dark shaded histogram shows the stars formed in the burst 
mode (see text). 


of the principal star formation recipes used in contemporary 
semi-analytic models, and we refer the reader to this paper 
for a comparison of the different techniques. 


4 SAMPLE MODEL GALAXIES 

For this study, we initially make use of the 450 realisations 
of model elliptical galaxies (bulge-to-total > 0.6) described 
in Beasley et al. (2002). These model galaxies span a B- 
band magnitude range of -22.7 < Ms - 5 log h < 19.5, and 
cover a range of environments (halo circular velocity) and 
masses. For direct comparison with the NGC 5128 MDF, 
we have selected eight model galaxies with luminosities that 
bracket the NGC 5128 luminosity of Mg ~ -21.0, and look 
for galaxies residing in groups rather than clusters. Beyond 
these luminosity and environmental constraints, the galaxies 
were chosen arbitrarily. 

We list some of the relevant characteristics of NGC 5128 
and the eight selected models in Table 1. The integrated 
B — V colours of the model ellipticals are very similar to 
that of NGC 5128^ ((B — 17)0= 0.88; de Vaucouleurs et al. 

T It should be noted that due to the large spatial extent of 


1991). The model V — K colours are somewhat redder than 
those given in Pahre (1999), although these were determined 
with unmatched apertures. This is agreement in the B — V 
colours is significant, since the fiducial Cole et al. (2000) 
model has difficulty in reproducing the red colours of lumi¬ 
nous spheroids (Cole et al. 2000; Benson, Ellis, & Menanteau 
2002 ). 

The output metallicity histograms for the total stellar 
populations in eight model galaxies are shown in Figure 2. 
In each case, the distributions reflect the relative number 
of stars formed per metallicity bin, for our adopted yield, 
feedback efficiency and IMF (in this case Salpeter, see Cole 
et al. 2000 for a full list of fiducial model parameters). 

The most luminous model galaxy (#1) is a central su¬ 
pergiant in a rich cluster, while the two faintest (#7 and 
#8) are galaxies smaller than NGC 5128. These are shown 
deliberately to help display the range of model outcomes 
that are possible. The other five are modef galaxies that are 
roughly similar in size to NGC 5128 and drawn from roughly 
Centaurus-like group environments. 

The different contributions of the quiescent and burst 

NGC 5128 on the sky, and the presence of dust in the disc of 
this galaxy, the reddening correction is somewhat uncertain (Is¬ 
rael 1998). 
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Figure 3. Age distributions for the stars in the eight sample models from Figure 2. Here, the ages are plotted as look-back times with 
the first star formation occurring at ~ 13 Gyr. The light filled histograms represent the star formation in the quiescent mode, the dark 
filled histograms represent the contribution from the burst mode. 

modes of star formation are plainly seen in these histograms: 
the ongoing quiescent mode, which takes place within the 
PGDs before they merge into the larger proto-galaxy, gen¬ 
erate all the metal-poor stars in the distribution. It is 
broad and nearly featureless, and its range and maxima are 
very similar from one galaxy to the next in our luminosity 
range. By contrast, the burst portion which is predominantly 
metal-rich can differ widely between galaxies: each merger 
produces a 5-function of stars at the metallicity of the col¬ 
liding gas, with amplitude in proportion to the amount of 
incoming gas available. 

In Figure 3, we show the age distributions (look-back 
times) for the same eight sample models. Here, it is obvious 
that both star formation modes may continue to surprisingly 
late epochs (though not necessarily so; some of the galaxies 
end up their formation quite early, particularly in the dense 
environments that give rise to rich galaxy clusters). Later 
infall of gas (and subsequent bursts) is more prevalent in 
lower-density environments. 

The galaxy-to-galaxy differences in the burst mode be¬ 
come particularly obvious for smaller elliptical galaxies. A 
smaller one which happened to be built by just one or two 
moderately large mergers will end up with a relatively nar¬ 
row metal-rich stellar component, whilst a galaxy that was 
built by a longer sequence of many smaller bursts will have 
a much broader final MDF but will often have a rather sim¬ 
ilar mean stellar metallicity. On the other hand, the biggest 


ellipticals are almost always the end result of a large num¬ 
ber of individual bursts, and so the metal-rich portion of the 
MDF for the biggest ellipticals is usually broad and moder¬ 
ately smooth once all these bursts are added together. 


Since the quiescently-formed stellar component is fairly 
similar amongst all the ellipticals in the semi-analytic model, 
it is the merger-formed metal-rich end which displays most 
variation from one galaxy to another. In essence, the metal- 
rich end of the stellar MDF should become less broad, more 
fragmented and with larger galaxy-to-galaxy dispersion as 
we go to lower-luminosity galaxies. 


In Figures 4 and 5 we show the MDFs and age distribu¬ 
tions for the GCs in these same sample models. In all cases, 
we have adopted the parameters es = 0.002, en = 0.007, and 
Ztrunc = 5 as described previously. Thus for the blue GCs 
in each system, the age distribution is very sharply peaked, 
happening entirely between 13 Gyr and 12 Gyr look-back 
time. It is only for the red GCs that late formation times 
occur. The MDFs are generally bimodal in form, although 
the metal-rich portions have more complex substructure. 
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Figure 4. Sample metallicity distribution functions for the GCs in the eight model galaxies plotted in the previous figure. The light 
filled histograms represent the star formation in the quiescent mode, the dark filled histograms represent the contribution from the burst 
mode. 


5 CONFRONTING THEORY WITH 
OBSERVATION 

To pursue a direct comparison between the models and the 
NGC 5128 data, we superimpose the predicted MDFs of the 
eight model galaxies on these NGC 5128 data. 

First, however, we need to account for the fact that 
the full information present in the model MDF will in ob¬ 
servational terms be partly lost by the smearing effect of 
measurement uncertainty in [m/H]. The effect is to signif¬ 
icantly smooth out the bursty, metal-rich end of the MDF 
while doing little to change the broad and already-smooth 
metal-poor end. For the NGC 5128 data from Harris & Har¬ 
ris (2002), the mean uncertainty is oy m / H ] = 0.06 dex.t 

We show the comparisons between the eight models the 
and NGC 5128 data in Figure 6. It is apparent that in their 
general characteristics, the agreement between the MDFs 
of the models and data is reasonably close. The shape and 
extent of the low-metallicity tails are generally similar, with 
small fractions of stars extending down to [m/H] ~ -2.6. 

t The measurement uncertainty in heavy-element abundance Z. 
as discussed by Harris &; Harris (2002), is equivalent to 0.06 dex 
scatter in [m/H] for the upper giant branch stars used in their 
study. Since the grid of giant-branch models has a track spacing 
of ~ 0.1 dex, this uncertainty corresponds to about half a grid 
step. 


There is no indication of an over-production of metal- 
poor stars ([m/H] < -0.7) in the model - ’the G-dwarf prob¬ 
lem’ - which is evident when a closed-box model is fit to 
these data (e.g. Harris & Harris 2002). In contrast to the 
closed-box ‘simple model’, where all heavy elements remain 
within the galaxy, the yield, y, in the semi-analytic model 
is a function of the level of metal ejection into the hot gas 
phase (e) and the feedback efficiency, and becomes an effec¬ 
tive yield, y eB : 


( 1 - 
1-R-/3 


(5) 


where R is the fraction of mass recycled by stars. In this way, 
material processed by stars is mixed into the halo gas reser¬ 
voir out to large radii. Therefore, large quantities of metal- 
poor stars are not produced in the semi-analytic model be¬ 
cause there is continuous infall of gas onto the star forming 
PGDs from this cooling reservoir (Kauffmann 1996). 

The metal-rich ends of the model MDFs in Figure 6 also 
show strong similarities to the NGC 5128 MDF. Both ex¬ 
hibit significant increases of stars towards high metallicity, 
peaking at or near solar abundances. However, it is at this 
metal-rich end that the differences are also most apparent. 
The model MDFs often extend to higher metallicities than 
is the case for the observed MDF. This yields mean metallic¬ 
ities between 0.1-0.2 dex more metal-rich than these data, 
depending upon the model galaxy selected. In some cases, 
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Figure 5. Age distributions for the GCs in the eight sample models, plotted as in the previous figure. The light filled histograms represent 
the star formation in the quiescent mode, the dark filled histograms represent the contribution from the burst mode. 


this excess is manifest as an excess of stars at Z > Zq, 
although for models #1, #5 and #8 in particular it also 
appears that the effective yield is possibly too high, since 
they over-predict the number of stars at ~ O.IZq. Further 
inspection also reveals that the excess of stars at the metal- 
rich ends of the models result in MDFs which are ~ 0.1 dex 
broader than these data. This is most evident for galaxy #1, 
the luminous gE galaxy in the centre of a rich cluster. 

In assessing these comparisons, it must be born in mind 
that the observed MDF almost certainly underestimates the 
numbers of high-metallicity stars within NGC 5128 in total 
for two reasons: (a) as is discussed in Harris & Harris (2002), 
the number of red giant stars with Z > Zq found in the 8 
kpc field is probably underestimated because of photometric 
incompleteness for these extremely red stars; and (b) the 
data samples only one location projected on the outer bulge. 

Real elliptical galaxies have radial abundance gradients 
(e.g. Davies, Sadler, & Peletier 1993; Gorgas et al. 1997), and 
so the inner bulge of NGC 5128 may contain more metal-rich 
stars. Both of these effects would raise the upper envelope of 
the observed MDF and increase the metallicity at the MDF 
peak. By contrast, the model galaxies represent a globally 
averaged population because they contain no spatial infor¬ 
mation. The average abundance gradient of giant ellipticals 
is ~ -0.2 ± 0.1 dex (e.g. Davies, Sadler, & Peletier 1993). 
This corresponds to a reduction of 40% over a factor of 10 
in radius, or a change in metallicity of 0.1 0.3 dex in going 

from 3 to 30 kpc. However, it remains to be seen whether 


this gradient will be manifest in a global metallicity shift in 
the observed distribution, or result in a change in the overall 
shape of the MDF. 

Bearing in mind the dearth of observed stars at Z > Zq, 
we perform Kolmogorov-Smirnov tests on the eight model 
galaxies to quantify their goodness-of-fit with these data. 
The resulting probabilities that the models and data are 
drawn from the same distribution, p(KS), are listed in Ta¬ 
ble 1. In terms of the model, as shown in Figure 2 the burst- 
formed stellar component displays a wide range of morpholo¬ 
gies, with some model realisations providing better fits than 
others. The best fit is achieved for galaxy #4, the elliptical 
with Mb - 5 log /i=-21.13, since the metal-rich end of this 
particular model is fairly narrow compared to other model 
galaxies in the sample. 

We compare a selected model galaxy (in this case #2) 
and its GC system with the NGC 5128 data in Figure 7. This 
galaxy has been chosen somewhat arbitrarily, although its 
luminosity and environment is not too dissimilar to that of 
NGC 5128. The excess of metal-rich stars is particularly ev¬ 
ident in this figure. For the GC comparison, we have shifted 
the Held et al. (2002) data by + 0.3 dex in order to ap¬ 
proximate a conversion between [Fe/H] and [m/H]. Both 
the model and data are clearly bimodal; the model peaks 
at [m/H] ~ -0.95 and -0.05, whilst the Held et al. dataset 
has peaks at [Fe/H] -1.2 and -0.3, corresponding to roughly 
[m/H] ~ -0.90 and 0. The metal-rich end of the GC MDF is 
coincident with that of the model stellar MDF because both 
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Figure 6. The globally averaged model MDFs (solid lines) for the eight selected galaxies compared to the NGC 5218 data (shaded 
histograms). The models have been normalised to contain the same number of stars as these data, and have been smoothed with an error 
corresponding to ~ 0.06 dex. 


have the same origin in merger-induced star formation. As 
discussed in Beasley et al. (2002), a good fit is achieved at 
the metal-poor end only if the GC formation in PGDs is 
halted early (z ~5). 


6 THE EFFECT OF STAR BURSTS ON THE 

MDF 

Since the primary purpose of this paper is to compare the 
predictions of the semi-analytic model with the observed 
MDF, our previous comparisons have been performed with 
the fiducial Cole et al. (2000) model. However, to test the ro¬ 
bustness of these predictions, we now change the relative im¬ 
portance of merger-induced star formation (the burst mode) 
in the model. The following is not meant to be an exhaus¬ 
tive exploration of the parameter space, but rather is a test 
of the sensitivity of the model MDF output to the burst 
assumptions. 

Two parameters which have a direct influence on the 
burst mode are / e iii P and /burst- As mentioned in Section 3, 
/eiiip represents a mass ratio between the central and satel¬ 
lite galaxy, which determines whether a merger is classed 
as ’major’ or ’minor’. Major mergers result in a change in 
galaxy morphology (i.e. discs are destroyed, stars are redis¬ 
tributed into the bulge) and in general any cold gas is con¬ 
verted in its entirety into stars. Minor mergers result in a 


mass accretion, but do not give rise to disc destruction, or 
generally, star formation. Any cold gas belonging to the in¬ 
coming satellite is added to the disc of the central galaxy. 

The fiducial model of Cole et al. (2000) sets / e iii P =0.3 
(a 30% central-satellite mass ratio), which is the lower limit 
of the range suggested by numerical simulations of mergers 
(Walker, Mihos, & Hernquist 1996; Barnes 1998). Increasing 
/eiiip has the principal result of changing the morphological 
mix of galaxies in the model; a higher mass ratio results in 
fewer major mergers and more ’discy’ systems (e.g. Benson, 
Ellis, & Menanteau 2002). The value of / e iii p does not have 
a significant effect on the time (or metallicity) of the last 
merger and its star formation, since the model ellipticals 
are generally selected on the basis that they have under¬ 
gone a merger recent enough to destroy any significant disc 
present. However, increasing / e m p does decrease the burst- 
formed component of the MDFs in the model, which is also 
governed by /burst- 

The parameter /burst sets the threshold for a merger to 
cause a burst of star formation, and this is set to equal / e m p 
in the fiducial model. In this way, a burst of star formation 
implicitly occurs whenever there is a major merger with cold 
gas present. We can test the effect on the mean galaxian 
MDFs by changing /burst to a smaller value, thus lowering 
the mass threshold required for a burst and increasing the 
mean number of bursts per galaxy. Setting /burst < /eiiip 
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Figure 7. Upper left panel: Comparison between a semi-analytic model stellar metallicity distribution function (solid line), and the 
observed MDF of old stars in NGC 5128 (shaded histogram). The model MDF has been smoothed by the measurement uncertainty of 
0.06 dex as discussed in the text. Upper right panel: Cumulative KS distributions of the model stars (dotted lines) and data. Lower left 
panel: Comparison between a model GC MDF (solid line) and the observed NGC 5128 GCs (Held et al. 2002; shaded histogram). The 
same smoothing as stated above has been applied to the model MDF. Lower right panel: Cumulative KS distributions of the model GCs 
(dotted lines) and data. 


allows impinging satellites to produce a burst of star forma¬ 
tion without disrupting any disc present. 

We have therefore re-simulated ~ 500 early-type galax¬ 
ies (B/T ratio > 0.6) in approximately the same luminosity 
range as the previous realisations (-22.7 < Mb - 5 log h 
< -19.5). For this purpose, we use the updated version of 
the semi-analytic code described in Benson et al. (2002). To 
facilitate a direct comparison with our previous model re¬ 
alisations, we have turned off the photoionizing background 
contribution (Benson et al. 2002) . The updated code pro¬ 
duces essentially identical integrated properties to those of 
our previous realisations for luminous ellipticals at 2 = 0. 

We compare the resulting MDFs with /b ur st=0.3 (left 
panels) and /burst=0.1 (right panels) in Figure 8. As shown 
in Section 5, whilst the form of the fiducial model MDF 
is very similar to these NGC 5128 data, the model is ~ 
0.1 dex offset towards the metal-rich end. The linear plots 
(Z/Zq) clearly show the excess of higher-metallicity stars. 
The model realisations with /b ur st=0.1 (i.e. have on average 
more bursts of star formation) are slightly more metal rich 
again since the stars formed in bursts which contribute to 
the model MDF all have around solar metallicities. However, 
the differences are not significant; within the scatter the two 
sets of model realisations are identical. 


In terms of the GCs, the model with /burst=0.1 produces 
proportionally more metal-rich GCs than the fiducial model 
used in Beasley et al. (2002). From comparisons with the 
well-studied elliptical NGC 4472, Beasley et al. (2002) found 
that the formation efficiency (in the absence of dynamical 
destruction) of the metal-rich GCs, £r=0.007. By apply¬ 
ing the normalisation to NGC 4472 described in Beasley et 
al. (2002), we find that £r=0.004 for /burst=0.1. Thus /burst 
effects the efficiency of the metal-rich GC formation, cr. 


7 AN AGE-METALLICITY RELATION FOR 
NGC 5128? 

We show the age-metallicity relation (AMR) of our model 
galaxy and its GCs in Figure 9. The dashed line in the 
figure indicates the average metallicity of stars at a given 
mass-weighted age in the galaxy. At 2 =0, this galaxy has 
a mass-weighted mean age of 9.5 Gyr, and a mass-weighted 
metallicity of [m/H]= -0.1. Prior to ~ 10 Gyr, the rela¬ 
tion for the galaxy stars is steep, and rapidly flattens off at 
later epochs (e.g. Kauffmann 1996). This is broadly consis¬ 
tent with the behaviour of the solar neighbourhood AMR 
(Edvardsson et al. 1993; Garnett & Kobulnicky 2000), and 
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Figure 8. Upper left panel: The mean model MDF produced by ~ 450 realisations of the fiducial model (with /burst=0.3) compared to 
the NGC 5128 data (shaded histogram). Error bars represent the 1-cr scatter on the model distributions. Upper right panel: The mean 
model MDF produced by ~ 500 realisations of the semi-analytic model with /'| )l]rst =0.1 (a minimum 10% central to satellite mass ratio 
is required to induce a star burst) compared to the NGC 5128 data (shaded histogram). Lower left panel: MDFs plotted on a linear 
scale. The shaded histogram represents the NGC 5128 data, the solid line shows the mean model distribution. Lower right panel: Linear 
version of the upper right plot. 


that of the SMC star clusters (Da Costa & Hatzidimitriou 
1998). 

The metal-poor GCs do not trace the enrichment his¬ 
tory of the galaxy, but it should be noted that this feature 
of the blue clusters is a direct consequence of the model 
assumptions (the formation epoch of these clusters is delib¬ 
erately truncated at z =5). However, they clearly are able 
to enrich substantially (up to [m/H]= 1.0) prior to our 

truncation redshift. This is a result of the rapid enrichment 
undergone by ’dwarf’ size haloes in the model. 

On the other hand, as metal-rich GCs are created in 
each major merger undergone by the galaxy, each burst of 
star formation produces GCs with ages and metallicities 
corresponding to the epoch of the merger. In this scheme, 
the metal-rich GCs are discrete tracers of the dissipational 
merging undergone by the galaxy. This is in contrast to the 
luminosity-weighted properties of the host galaxy since, not 
only will age-fading make any burst-formed stellar popu¬ 
lation virtually undetectable after a few Gyr (e.g. Bruzual 
& Chariot 2001), but also the rapidly decreasing cold gas 
fraction in PGDs towards later epochs leads to increasingly 
smaller burst-formed populations (e.g. see figure 2 in Beasley 
et al. 2002). 

It is also worth noting in Figure 9, that a fraction of the 


metal-rich GCs are slightly more metal-rich than the global 
properties of the galaxy itself. This occurs because, whilst 
the GCs are formed from enriched, well-mixed gas of a sin¬ 
gle metallicity, the galaxian stellar population comprises of 
a mixture of metal-rich and metal-poor stars (e.g. Figure 7). 
Since merging is expected to produce new stars preferen¬ 
tially in the central regions of the merger product (e.g. Mihos 
& Hernquist 1996), the central kpc of the galaxy may be 
expected to have properties most similar to the mean prop¬ 
erties of the metal-rich GCs. 

Note that the extent of the blue GC AMR is sensitive 
to the truncation redshift we adopt, truncating their forma¬ 
tion at later epochs results in more metal-rich GCs. The red 
GC AMR is a direct result of the merger histories in the 
hierarchical model. 


8 DISCUSSION AND CONCLUSIONS 

Semi-analytic models have had some success in reproducing 
the integrated properties of galaxies over a wide distribu¬ 
tion in redshift (e.g. Kauffmann, White, & Guiderdoni 1993; 
Kauffmann 1996; Baugh, Cole, & Frenk 1996; Baugh et al. 
1998). In this paper, we have pursued a rather more spe- 
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Figure 9. The predicted age-metallicity relation (AMR) for our 
NGC 5128-like model galaxy and its GCs. The dashed line indi¬ 
cates the AMR for the galaxy, which at z =0 possesses a mass- 
weighted mean age of 9.5 Gyr, and a metallicity of 0.8 Zq. Filled 
circles are the (binned) metal-poor GCs, which have had their 
formation truncated at z =5. These have very similar ages (~ 
12.4 Gyr), but are able to substantially enrich at early times. 
The metal-rich GCs are shown as filled triangles, and mark the 
epoch of each dissipative merger undergone by the galaxy. Their 
formation proceeds to late epochs in this galaxy. 



cific comparison; that between the observed and predicted 
metallicity distribution functions of the old stars and GCs 
of a single elliptical galaxy. 

Clearly some discrepancies arose in the comparison; as 
noted above, the most obvious one is that the metal-rich stel¬ 
lar component of the model is broader and more metal-rich 
than the observed MDF. However, it is important to remem¬ 
ber we have compared an ensemble of model realisations to 
data for a single galaxy. Whilst integrated colours and spec¬ 
tra give some reason to believe the MDFs of ellipticals may 
be similar, the MDF of the outer bulge of the Milky Way 
(Zoccali et al. 2002) shows some subtle differences from that 
of NGC 5128. Until more star-by-star MDFs are obtained 
for other ellipticals, the existence of significant variations (or 
lack thereof) in their MDFs has yet to be demonstrated. 

In terms of the model, we could have altered a number 
of parameters to better fit the galaxy MDF. For example, de¬ 
creasing the yield in the model by 0.1 dex would act to lower 
the global metallicity of the model MDF, better matching 
the metal-poor tail of the observed MDF (and GCs) shown 
in Figure 7. Since the expected theoretical yields for a given 
IMF are rather uncertain (e.g. Prantzos 2000) this would 
not be an unreasonable adjustment, although in so doing 
may weaken the successes of the model in reproducing other 
galaxy properties. 

In any event, obtaining an exact match to the data was 
not the purpose of this study. Rather, our purpose was to 


see whether the star formation in the model leads to an 
adequate reproduction of the galaxy data, and whether in 
turn we may reconstruct a plausible star formation history 
for NGC 5128. In this respect, the agreement has been en¬ 
couraging. The vast majority of the star formation in the 
model occurs ’quiescently’, where the interplay between gas 
cooling, star formation and feedback yields old stellar pop¬ 
ulations which are chemically enriched not unlike those ob¬ 
served in the NGC 5128 bulge/halo. Merging of this enriched 
gas creates an additional metal-rich component, including 
metal-rich GCs (e.g. Beasley et al. 2002). 

Considering the relative simplicity of the star forma¬ 
tion prescription in the semi-analytic model, the agreement 
between the observed and model MDFs is remarkable. We 
argue that this is also true for the metal-rich GCs, which fall 
naturally out of a hierarchical merging model. Constraining 
the red GC formation efficiency through observation is a 
reasonable undertaking until we have a consistent theory 
for globular cluster formation (e.g. Ashman & Zepf 2001). 
Formation of the blue GCs in the model requires for them 
to be disconnected from the star formation in their progen¬ 
itors after the truncation redshift Ztrunc- Understanding the 
formation of the blue GCs will require an understanding of 
these high-redshift progenitors (e.g. Bromm & Clarke 2002; 
Carlberg 2002). 

To conclude, we find that the metallicity distribution 
functions of stars and globular clusters in the nearby el¬ 
liptical galaxy NGC 5128, currently the only such galaxy to 
have its MDF derived on a star-by-star basis (Harris, Harris, 
& Poole 1999; Harris & Harris 2000; 2002), are well repro¬ 
duced by a semi-analytic model of galaxy formation (Cole 
et al. 2000; Beasley et al. 2002). Since the imprint of rela¬ 
tively old mergers may be hidden in stellar populations of 
galaxies, GCs may provide one of the few probes of this star 
formation. A large spatially unbiased sample of good quality 
spectra for NGC 5128 GCs will provide an important test 
of this scheme. 
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